arXiv: 1503.00379v2 [cond-mat.str-el] 31Jul2015 


Renormalization group analysis of the pair-density-wave and charge order within the 
fermionic hot-spot model for cuprate superconductors 

Hermann Freire^0 Vanuildo S. de Carvalho^’^, and Catherine Pepin^ 

^ Instituto de Ftsica, Universidade Federal de Goids, 74-001-970, Goiania-GO, Brazil and 
^ IPhT, L’Orme des Merisiers, GEA-Saclay, 91191 Gif-sur-Yvette, France 
(Dated: August 4, 2015) 

In light of the new experimental and theoretical important developments in high-Tc superconduc¬ 
tivity, we revisit the fermionic hot-spot model relevant to the phenomenology of the cuprates. We 
extend previous results by means of a complete two-loop order renormalization group (RG) frame¬ 
work. Here, we explicitly study the effect of the charge-density-wave (CDW) order parameter with 
a d-wave form factor with the experimentally observed modulation (±(3o,0) and (0, ±Qo) at the 
infrared-stable nontrivial hxed point obtained previously for this model. Additionally, we proceed 
to investigate also the so-called pair-density-wave (PDW) order that was recently proposed in the 
literature as a possible candidate for the “hidden” order to describe the pseudogap phase observed 
in underdoped cuprates. We conhrm that although the above two ordering tendencies are also found 
to be nearly degenerate both at one-loop and two-loop RG orders and linked by an emergent SU (2) 
pseudospin symmetry, they turn out to be subleading for weaker couplings in the present model to 
antiferromagnetism, d-wave bond-density wave (BDW) order with modulation along Brillouin zone 
diagonals {±Qo,±Qo), and d-wave singlet superconductivity (SSG). However, as we increase the 
strength of the initial coupling towards moderate values, we do capture a tendency for the entan¬ 
gled PDW/CDW order to become leading compared to BDW/SSC in the model, which suggests 
that the former composite order might be indeed a viable concept to describe some cuprate super¬ 
conductors at high temperatures in the underdoped regime, as has been recently alluded to by many 
authors in the literature. 
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I. INTRODUCTION 

The underlying nature of the pseudogap phase ob¬ 
served in the underdoped cuprates continues to be one 
of the most enigmatic and profound problems in con¬ 
densed matter theory. Despite this statement, a lot has 
been learned about this phase both theoretically and 
experimentally in the last years. On the experimen¬ 
tal side, groundbreaking measurements performed on the 
non-Lanthanum-based materials have increased dramati¬ 
cally our present knowledge about these compounds. Nu¬ 
clear magnetic resonance experimentsii^, resonant^id and 
hard^ x-ray scattering, ultrasound measurements^ and 
scanning tunneling microscopy^i^ confirmed the emer¬ 
gence of a charge-density-wave order (CDW) at small 
hole doping at wavevectors (±(5o,0) and (0, ±Qo) ori¬ 
ented along the principal axes of the Cu 02 planes with 
a predominant d-wave form factor— dS. This CDW has 
the form of a checkerboard and is incommensurate with 
the lattice, since it connects approximately the “tips” 
of the Fermi arc o^^d^ . This gave an initial impetus for 
the interpretation that such a CDW could be the miss¬ 
ing “hidden” order in the pseudogap phase that might 
be ultimately responsible for the reconstruction of the 
Fermi surface observed, e.g., in quantum oscillations ex¬ 
periments in YBCO at high magnetic field a^^d^ . How¬ 
ever, this latter interpretation, at least at face value, has 
been recently challenged by some experiments due to the 
fact that the transition line of this CDW phase is essen¬ 
tially dome-shaped and, in general, is not coincident with 
the pseudogap line T* observed in the phase diagram of 


these compounds. 

Since the “tips” of the Fermi arcs in the pseudogap 
phase are close to the so-called hot spots (i.e., the points 
in momentum space in which the antiferromagnetic zone 
boundary intersects a putative underlying Fermi surface 
of these compounds), this could mean that a possible 
low-energy effective model to describe these materials 
might be a spin-fermion model, or more generally, a hot¬ 
spot modeli^^— . In this respect, it has been famously 
demonstrated in Ref— ^ that, if the energy dispersion is 
linearized, an emergent SU(2) pseudospin symmetry that 
rotates the d-wave superconducting order parameter onto 
a “d-wave bond-density wave” (BDW) order— with mod¬ 
ulation along Brillouin zone diagonals (iQojiQo) ex¬ 
ists in the model and this property effectively produces 
a composite order parameter with both bond order and 
preformed pairs at high temperatures^^. However, a gen¬ 
erally acknowledged drawback of this approach is related 
to the fact that the leading charge order indeed always 
appears at the wavevector (±(5o,±Qo)i which, so far, 
has never been observed experimentally. Many alterna¬ 
tive scenarios22r.2£ that aim to resolve this discrepancy 
with experiments have also been proposed in the litera¬ 
ture, but no consensus has yet emerged in the community 
with regard to what is the mechanism responsible for the 
pseudogap phase in the cuprates. 

For this reason, many researchers continue trying to 
uncover what might be the “hidden” order that is fun¬ 
damentally responsible for the emergence of the pseu¬ 
dogap phase seen in the cuprates. Several good candi¬ 
dates on the table (which are not mutually exclusive) in- 
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elude pair-density-wave (PDW) order— the aforemen¬ 
tioned d-wave CDW order at wavevectors (iQoi 0) and 
(0, ±(5o) (which may also break additional discrete sym¬ 
metries such as the (74 lattice rotational symmetry down 
to the C 2 , and time-reversal and parity symmetrie a^^i^^ ), 
Mottness^i^, orbital loop current orderi^^^^l^ fractional¬ 
ized Fermi liquid^l^ among others^^. For simplicity, we 
shall concentrate in this work only on the first two candi¬ 
date orders (i.e., PDW and CDW) and leave the analysis 
of other possible orders that might be present in the un¬ 
derdoped cuprates for a future study. 

The PDW order refers essentially to a superconducting 
order with a finite Cooper-pair center of mass momen¬ 
tum, similar in this respect to a Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state^>^. It has been recently 
proposed that this order may account for the anoma¬ 
lous quasiparticle excitations seen by angle-resolved pho¬ 
toemission (ARPES) experiments in both LBCO and 
YBCO^. The PDW order could also potentially emerge 
from a doped Mott insulator scenario as a strong¬ 
coupling instability due to an effect similar to “Ampere’s 
law”—. In this respect, in Ref— it was argued that 
the pseudogap phase in underdoped cuprates is more 
appropriately described by a PDW instead of simply 
a CDW. Additionally, it has also been pointed out in 
the literature that the PDW order might give rise to 
a secondary non-superconducting order parameter that 
breaks both time-reversal and parity symmetries, but 
preserves their product—. The breaking of time-reversal 
and parity symmetries in the pseudogap regime would 
be consistent, e.g., with the experimental signatures ob¬ 
tained both via polarized neutron scattering^ and Kerr- 
rotation experiment s"*^^ 7^ . Therefore, it becomes clearly 
important to analyze, in a more detailed way, the role of 
the PDW order parameter in the context of low-energy 
effective models that might describe the essence of the 
physics of the underdoped cuprates. 

In light of the many recent experimental and theoreti¬ 
cal developments described above, we build upon a work 
by two of us by revisiting here the so-called fermionic 
hot spot model defined in Ref— and analyzed within a 
two-loop renormalization group (RG) framework. The 
fermionic hot-spot model can be seen as descendant of 
the Abanov-Chubukov spin-fermion modeU^ that, most 
importantly, includes here all relevant interactions be¬ 
tween the fermions and, for this reason, enables one to 
analyze all possible instabilities of the model from weak 
to moderate couplings on equal footing. We mention that 
the fermionic hot spot model has also been recently stud¬ 
ied by Whitsitt and Sachdev^ within a one-loop RG ap¬ 
proach and their results, up to this order, agree with 
ours. In the present paper, we explicitly study the ef¬ 
fect of the charge-density-wave (GDW) order parameter 
with a d-wave form factor with the experimentally ob¬ 
served wavevectors {±Qoi 0) and (0, ±(5o) in the vicinity 
of the infrared-stable nontrivial fixed point obtained at 
two loops^ that, as a result, naturally implies a new 
quantum critical universality class for the present model. 



FIG. 1. (Color online) The present 2D fermionic model con¬ 
sisting of eight hot spots (with bare Fermi velocities vf and 
angle 6) on the Fermi surface that are directly connected 
by the commensurate spin density wave (SDW) ordering 
wavevector Q = (7r,7r). The incommensurate wavevectors 

Q 3 : ((5o,0), (^y (O 5 Q 0 ), Q (( 5 o,( 5 o), Q 3 ; and Qy are 

also displayed. 


Moreover, we proceed to investigate in this scenario also 
the role of the PDW order, which was recently proposed 
in the literature as a viable candidate for the “hidden” 
order to describe the pseudogap phase observed in the 
underdoped cuprates. Interestingly, we confirm that the 
above two orders turn out to be also nearly degenerate 
both at one-loop and two-loop RG orders and clearly 
linked by an emergent SU (2) pseudospin symmetry at 
low energies. Despite this fact, they continue to be sub¬ 
leading for weaker couplings in the present model to anti¬ 
ferromagnetism, d-wave bond-density wave (BDW) order 
with modulation given by {±Qo, ^Qo) and d-wave singlet 
superconductivity (SSG). We then finish by discussing 
how the PDW/GDW order is likely to become leading 
compared to BDW/SSC order in this model, which would 
potentially agree with the experimental situation. 

This paper is organized as follows. In Section II, we 
define the so-called fermionic hot spot model that, as we 
will explain, might be relevant to describe some aspects 
of the phenomenology of the underdoped cuprates. In 
Section III, we briefly introduce the field-theoretical RG 
approach and show how to implement this method up 
to two loops in the present model. In this part, we will 
concentrate specihcally on discussing the role played by 
the PDW and GDW orders (both with a predominant 
d-wave factor) at the experimentally observed wavevec¬ 
tors (±Qo, 0) and (0, ±Qo) that emerge in the low-energy 
limit of the model. The relevant RG flow equations are 
then derived and solved in Section IV. Finally, we present 
our conclusions in Section V. 
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FIG. 2. The relevant couplings in the present fermionic hot-spot model. For the gis, gst and gsv couplings, we show two 
possible scattering processes (one by a solid arrow and the other by a dashed arrow), which are always equal in our two-loop 
RG results. 


II. THE FERMIONIC HOT-SPOT MODEL 

To begin with, we consider a single-band fermionic 
model whose noninteracting part is given by the fol¬ 
lowing energy dispersion = —2t(cos(fca,) -|- cos(fcy)) — 
At' cos{kx) cos(fcy) — /r, which is given in terms of units of 
the lattice constant, the parameters t and t' are, respec¬ 
tively, the nearest neighbor and next-nearest neighbor 
hoppings and /i is the chemical potential. For the spe¬ 
cific case of the cuprates, one may assume for simplicity 
the choice of parameter t' = —0.3t, which, at low hole 
doping, results in the Fermi surface (FS) displayed in 
Fig. 1. In the interval 0 < |p| < 4|t'|, the correspond¬ 
ing FS intersects the antiferromagnetic zone boundary at 
eight points (the aforementioned hot spots). These hot 
spots are shown in Fig. 1 and the excitations residing 
near the locus of these points in fact generate the most 
singular contributions in the model at long distances and 
long timescales. Besides, one can naturally rotate the 
momentum axes by an angle of 7r/4 by defining new axes 
{kx,ky), where the rotated momenta kx and ky are de¬ 
picted in Fig. 1. 

Since we will focus on the universal properties of this 
model, we may linearize the energy dispersion around 
the hot spots as ~ vp.k with k being the mo¬ 
mentum measured relative to the hot spots, the Fermi 
velocity 'VF={vx,Vy) with the modulus of the compo¬ 
nents given by Vx = |V(fc^_fc^)^k|k=kFl cos6> and Vy = 
l^(fcx.fc„)^k|k=kFl sin0, with 9 being the angle depicted 
in Fig. 1. Here, we specialize to the value /i « —0.77t 
(i.e., a low hole-doped regime). It is interesting to point 
out, however, that our results are fairly insensitive to this 
choice of the chemical potential, as long as 0 < \^\ < 4|t'|. 
The components of the momentum kx and ky are re¬ 
stricted to the interval [—fee, fee], where k^ stands for a 


sharp ultraviolet (UV) momentum cutoff in our theory. 
This implies also a high-energy cutoff, which is given by 
Aq = 2vpkc that we choose to be of the same order of 
the full bandwidth of the problem, i.e., Aq ^ 8t. 

The fermionic hot spot model at zero tempera¬ 
ture is described by the partition function Z = 
J'D[tj},ip]exp{i J^^dt Lii[tj;,tp]) with the interacting 
renormalized Lagrangian Lji given by 

Lr = Y,Z^ (k) [^^t - Z^v^kx - Zyy^ky] V^^(k) 

k,{7 

“X! X! ^|5*^V'a(k4)^/’f'(k3)V’a'(k2)V'^(ki), 

i ki,k2,k3 
(7,g' 

( 1 ) 

where = ki -|- k 2 — ka and the volume V has 
been set equal to unity. All bare quantities (indi¬ 
cated by the superscript B) are related to the renor¬ 
malized quantities (indicated by the superscript R) by 
the following standard expressions: (k) = z]p'^'tp^{]<i), 

V’f(k) = ^yVf(k), uf = Z^v^ and 

where Z^ is the field renormalization factor and 

and Zy are the Fermi velocity renormalization factors. 

The renormalized Grassmann fields V’o-(k) and ^^(k) 
are associated, respectively, to the creation and an¬ 
nihilation operators of the low-energy excitations re¬ 
siding in the vicinity of the hot spots with momen¬ 
tum k and spin projection a. The index i refers to 
the many possible interaction processes in the model 
that generate singularities within perturbation theory. 
Here we follow roughly the notation of Ref4^, i.e., i = 
1, 2,3, Ic, 2c, la:, 2x, Is, Ir, 3a:, 3p, 3<, 3u, 16, 2/ (for details 
of the couplings taken into account, see Fig. 2). In 
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this way, our approach will mirror the so-called “g- 
ology” notation^, adapted to our 2D problem at hand. 
Additionally, we define dimensionless renormalized cou¬ 
plings (denoted by gf) in the following way: gf = 
[g(^ + Sgf^], where Nq = kc/ir'^Vp is the den¬ 
sity of states at the Fermi level. In all above expressions, 
we set conventionally Z^p = l + 6Ztjj, Z^ = Z^\l-\-5Z^), 
Zy = Z^^{1 + SZy), where SZ^, SZ^, SZy and Sgf^ are 
the counterterms that will be calculated order by order 
within the renormalized perturbation theory^. 

Since the number of couplings in the present model 
is large, the possibilities for choosing their initial condi¬ 
tions are also considerable. Therefore, our choice here 
will be motivated by relevant physical microscopic mod¬ 
els. We shall restrict our analysis in the present paper to 
the paradigmatic 2D Hubbard model, whose initial con¬ 
ditions for the couplings are given by g^ = kJJ/{-K'^Vp) 
for all i, where f/ > 0 is the local on-site repulsive inter¬ 
action strength. 


III. FIELD-THEORETICAL RG 


In this section, we describe the two-loop RG approach 
that we shall apply in the rest of this work. Our dis¬ 
cussion here will be relatively concise, since the field- 
theoretical RG methodology was already explained by 
two of us in great detail in the context of the present 
model elsewhere^i^. Within perturbation theory, if we 
compute many one-loop quantities, such as, e.g., particle- 
particle and particle-hole polarization bubbles for several 
choices of incoming external momenta q and also their 
corresponding two-loop-order corrections, we obtain sev¬ 
eral logarithmic divergences of the type ln(Ao/A) as we 
probe the system towards the low-energy limit A —>■ 0 (see 
also Refs4^^— ). As we have already explained above, 
we circumvent this problem by defining counterterms 
for all the dimensionless couplings, the fermionic fields 
and the Fermi velocity of the model that, by construc¬ 
tion, effectively regulate all divergences order by order 
within the renormalized perturbation theory. The result¬ 
ing coupled-differential RG equations obtained by this 
method are thus given in Appendix A. 

To exemplify the field-theoretical RG method, if we 
calculate the self-energy of the present model (e.g., for 
kx > 0 and ky > 0) at two-loop order within the renor¬ 
malized perturbation theory, we obtain the following non- 
analytic contribution 


(2 loops) / 


R (k, W) = - ^^2(Ir^2 + IVyVyky] 


E 

X In 


An 


max{a;, uj^lkl} 


- [5Z^ uj - SZ^v^kx - SZy^ky] , 

( 2 ) 


superscript R in the above renormalized dimensionless 
couplings of the model from here on. As a result, we 
obtain the following RG flow equation up to two loops 
associated with the renormalized Fermi velocity = 
iVx,Vy), i-e., 


A 


dRR 

Ha 


Rr 


4(1 + 4 ) 




( 3 ) 


where kr = (Vy/v^) = tan0/j is the ratio between the 
components of the Fermi velocity at the hot spots and 
Hyiigf'}) = Ivy + lijj- From Eq. ([S]), one obtains that 
the only infrared (IR) stable fixed point corresponds to 
Kp ^ 0 (or 0p 0), which refers to a perfect nest¬ 
ing condition of the hot spots connected by the antifer¬ 
romagnetic (AF) ordering wavevector Q = (tt, tt). Al¬ 
though this only happens, strictly speaking, in the limit 
of A —0, one can check numerically that for moderate 
couplings and not too low RG scales, k is already rea¬ 
sonably suppressed and the resulting renormalized FS 
exhibits an approximately nesting condition connected 
by (7r,7r) (see, e.g.. Fig. 3). This fact will be clearly fa¬ 
vorable to antiferromagnetic fluctuations that will play 
a major role in the present system as will become clear 
shortly. 

It is important to discuss at this point on the contri¬ 
bution of both particle-particle and particle-hole bubbles 
for some choices of wavevector q (to check our present 
definitions, see Fig. 1). If we calculate those quantities, 
we obtain the following expressions 


TT t r\\ -^0 , 4o 

IIp_fT(q = Q) = — In 


A / ’ 


npp(q= Q) = -^In ) , 


IIpp(q = Q) 


Nn 


■In 


An 


npp(q = Qa:) ~ A^oAln 


2cos0p \max{A, 2 p^A:c sin0p} 




( 4 ) 

( 5 ) 

)' 

( 6 ) 

( 7 ) 


npp(q = Qa:) ~ —NoXln 


An 


niax{A, ’ 


( 8 ) 


IIpp(q = Q^) ~ iVoA'ln 


An 




(9) 


IIpp(q = Q^) Ri —A^oA'ln 


An 


^^ax{A, 


( 10 ) 


where 7V, = = {glc+glc+9lx+92x+9lp+9lx-9icg2c- 

5ix52x-53p53a,) and7„^ = {gfx+92x+93p+93x-9ic-92c+ 
gic 92 c-gixg 2 x- 93 pg 3 x)- For simplicity, we shall omit the 


with A being the RG scale which could be interpreted 
(to logarithmic accuracy) in terms of an effective tem¬ 
perature scale. Besides, the prefactors A and A' are 
















5 



FIG. 3. (Color online) A possible renormalized Fermi sur¬ 
face suggested by the present two-loop RG calculation of the 
fermionic hot-spot model. Note the resulting nesting condi¬ 
tion of the hot spots connected by Q = (tTjTt). 


given by A = l/(cos0fl; -|- sin^/j) and A' = l/(cos0i{ — 
sinOfi) for the present hot spot model. We mention 
here that even though the bubbles npjj(q = Q^.) and 
npp(q = Qa;) have initially different prefactors and dif¬ 
ferent low-energy cutoffs compared to npp(q = Q^) and 
npp(q = Q^), those prefactors and cutoffs become iden¬ 
tical for the renormalized Fermi surface (i.e., dp 0). 
In this way, we will refer only to one type of particle- 
particle and particle-hole bubbles [i.e., npp(q = Q^,) 
and npp(q = Qa,)] in what follows. All the above po¬ 
larization functions agree precisely with those derived in 
Ref4^. Those particle-hole and particle-particle bubbles 
at moderate energy scales satisfying A > 2tipfccSindp 
and A > Vpkc{l — tan dp)/A (which would correspond 
physically to a “high-temperature” regime) will approx¬ 
imately generate logarithmic divergences as a function 
of A that will allow us to use the field-theoretical RG 
approach to this model. A similar result also holds for 
the particle-particle and particle-hole bubbles at Qy = 
(OjQo) and Qy. For moderate interactions, it can be 
shown numerically that a nontrivial fixed point can be 
reached within such a high-energy regime. We empha¬ 
size here that such an approximation resonates with the 
idea that there is possibly a ‘high-temperature’ pseudo¬ 
gap (T < T*) in the cuprate superconductors, which 
plays the role of the ‘normal’ state in these compounds 
out of which many broken symmetry states (such as, e.g., 
a phase with either incommensurate charge order or pos¬ 
sibly also pair-density-wave order) emerge at lower tem¬ 
peratures (T <T** <T*). 

We first begin by reviewing previous results^ and then 
we present here new results concerning the present model. 
We integrate numerically the RG flow equations for the 
model [Eqs. (A1)-(A15)] using conventional fourth-order 
Runge-Kutta method. In one-loop RG order, the nu¬ 
merical solution of those equations shows that, despite 
the fact that there is no evidence of a non-trivial fixed 
point in the model at this order, nearly all dimension¬ 


less coupling constants diverge at the same critical RG 
step value 4 = ln(Ao/Ac), which depends only on the 
initial conditions for the couplings. This divergence is 
appropriately described by the following scaling ansatz 
gi{l) = Ci/(lc — l), with Ic oc (1/G) and the Ci are univer¬ 
sal constants that do not depend sensitively on the values 
of the initial conditions. For instance, for Hubbard-like 
initial conditions, their numerical values are found to be 


(Cl, C2, C3, Cic, C 2 c, Cix, C2X, Csx, Csp, Cst, C^v) ~ 
(0.07,0.23,0.39, -0.04, -0.02,0.13,0.19,0.19,0.11, 
0.12,0.25), (11) 

and Cis = Cir = Cib = C 2 / = 0. It is interesting to note 
that in the resulting effective theory there are indeed only 
four coupling constants that do not diverge {gis, gir, gib, 
and g 2 f) as the RG scale A is lowered towards Ac. In this 
way, they become comparatively much smaller than the 
other couplings at low energies. For this reason, one can 
go back to the original model [Eq. ([T|)] and neglect for 
simplicity all these irrelevant couplings to begin with. We 
shall perform this in the remainder of this work. 

Next, we move on to the two-loop RG case. An impor¬ 
tant point we wish to stress here is that it is crucial to 
implement the RG method for the present model at least 
at two loops or beyond, since as we have shown previously 
the first nonanalytic contribution to the self-energy (and, 
consequently, the renormalized single-particle Green’s 
function and related quantities^i^) only emerges at this 
order. In this respect, by solving numerically the com¬ 
plete two-loop RG equations [Eqs. (A5)-(A15)], we can 
show that one finds instead a “landscape” of many non¬ 
trivial fixed points in the model, in addition of course 
to the trivial (i.e., Eermi-liquid) one. Interestingly, only 
one nontrivial fixed point^ turns out to be stable in all 
directions in the IR regime, so that any trajectory in the 
coupling parameter space close to this fixed point will 
necessarily converge to it as the RG scale A 0. For 
this reason, we shall focus our attention from here on 
only on this single IR-stable nontrivial fixed point, since 
it naturally implies a new quantum critical universality 
class for the present problem. The numerical results for 
this nontrivial fixed point obtained at two-loop order are 
approximately given by 


5 92 5 93^ 9lc^ 920 9lx^ 92 x^ 93x 5 93p^ 93t 5 93v) ^ 

(0,1.68,1.84, -2.0, -1.0,1.92,1.92,1.92,0,0,0). (12) 


We will refer to the above nontrivial fixed point hence¬ 
forth as simply the hot-spot fixed point (HSFP). It is 
important to emphasize that Umklapp interactions in 
the present fermionic hot-spot model are fundamental for 
stabilizing the HSFP in the IR-limit. It may also be note¬ 
worthy the fact that some couplings, due to two-loop- 
order quantum fluctuations, now flow asymptotically to 
zero at low energies, whereas others flow instead to finite 
values. This interesting fact should potentially simplify 
the solution of the present model by means of other com¬ 
plementary methods. 
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IV. LINEAR RESPONSE THEORY 

To study the enhanced correlations and possible order¬ 
ing tendencies in the present model, we must add to Eq. 
(1) the following term 

Lext = q) + Cl/2)^B^{-k + q/2) 

k,a,/3 

+ ^BW(k,q)V'Ba(k + q/2)V'B/3(k-q/2) 
k,a,/3 

(13) 

where A^^sp(k,q) and (k, q) correspond to the 

bare response vertices for the superconducting (SC) and 
density-wave (DW) orders, respectively. This additional 
term will generate some new Feynman diagrams with 
three-legged vertices (see Fig. 4), which may also exhibit 
logarithmic divergences as a function of the RG scale 
A in the model. We can then conventionally define the 
renormalized response vertices with their corresponding 
counterterms as follows: A^^-(k, q) = Z^^[A^^-(k, q) -|- 

(5A^^-(k, q)] for i = SC and DW. As a result, we shall 
obtain straightforwardly the corresponding RG equations 
up to two loops. Since two of us have already derived 
elsewhere^ the RG equations for the order parameters 
associated with antiferromagnetism, d-wave supercon¬ 
ductivity and d-wave bond-density-wave order at the in¬ 
commensurate wavevector Q = {Qo,Qo), we shall not 
repeat them here. For this reason, we concentrate only 
on the PDW and d-wave charge order response vertices 
with modulation Q^, = (Qo, 0) [or Qy = (0, Qo)] in what 
follows. If we take into account the effect of the renor¬ 
malization of the FS, the corresponding RG equations 
are given by 



FIG. 4. (Color online) Feynman diagrams for (a) SC and 
(b) DW response functions of the present fermionic hot-spot 
model. Solid lines represent noninteracting fermionic single¬ 
particle Green’s function, while the wavy lines stand for the 
renormalized coupling constants. Depending on the channel 
analyzed, the vector q can be: Q = (tTjTt), Q = (Qo,Qo), 
Qx = (Qo,0), Qy = (0,Qo), QL and Q^. 

f A (j) _ A Wtl _ A Oft 

I ^PnW(Q^(y)) - ^E,SC(Q^(y)) ^Ii,SC(Q,(y)) 

] A O') _ AO)tt , A 0)44. 

where the points at the FS are denoted by j = 1,2 (see 
Fig. 5). To determine the symmetry of the order parame¬ 
ter, we must further symmetrize (or antisymmetrize) the 
response vertices with respect to the index j. Thus, we 
obtain 

{ A - _ a(2) 

- ^PDW{Qa:(y^) ^PDW{Qa:(y^)- 
J. (d-wave) _ A (1) _ A (2) 

^CDW(Ci^^y)) - 

.(s-wave) _ a(1) i A (2) 

^CDWiQ^l^y^) - ^CDW(Q^^y^) ^CW(Q„(j,))' 


A 




(.y )) 


dA 


= A 


gsv 


E * (2)aa 

^R,DW (Q^l^yf) 


- gstA 


(2)/3a 

R.DWiQ^^y-,) 


+ gA 


(l)«/3 

R,DW(Q^^y-,)’ 


(14) 


A 


dA 


(l)a/3 

R.SCjQ^^y^) 


dA 


= A 


gat A 


(2)a/3 


-I- »7 A 


(l)a/3 

R,SC(Cl^(^yf)’ 


gavA 


(2)/3a 

R,SC(Q^^y>i) 


(15) 


where the two-loop self-energy feedback contribution is 
given by the anomalous dimension rj = A(dlnZ^/dA) = 

{g1 + fli + gL + flic + flL + aL - flifl2 - flicfl2c - gixg2x - 
fl 3 pfl 3 a;+flip + flL + fl3/2)/4. It is interesting to point out 
here that the above RG equations are the same for both 
Qa; and Qy, such that the C 4 lattice rotational symme¬ 
try is always preserved in our theory. This explains our 
present notation Qx(y) in the response vertices. 

By antisymmetrizing these response vertices with re¬ 
spect to the spin indices, we obtain the PDW and CDW 
order parameters 


As a result, the two-loop RG flow equations for the re¬ 
sponse functions associated with a potential instability of 
the normal state towards a given ordered (i.e. symmetry- 
broken) phase then finally read 


' 'I ^pnw(Lixty))^ y 

A ^ A d—wave _ \ / o i \ A d—wave 

^dA^'^^^(‘4x(s)) “ '^^~^9av + gatj^cDWiQ^^y^) 

I _ \ d—wave /i-tA 

+ (17) 


A. One loop RG 


To get a qualitative idea of the results in the RG anal¬ 
ysis, it can be very useful here to discuss first the results 
of the present model in a one-loop approximation, before 
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FIG. 5. (Color online) The signs of the order parameters at 
all eight hot spots for a d-wave form factor. Instead, for an 
s-wave form factor, all signs become positive. 


moving on to the full two-loop RG case. For this reason, 
we shall perform this now. At one-loop RG order, we can 
calculate the PDW and GDW response functions close 
to the critical scale by simply substituting the ansatz 
gi{l) = Cij{lc — 1) for the couplings into Eqs. (fT6|) - (fT^ 
and setting the two-loop anomalous dimension contribu¬ 
tion to zero (i.e., rj = 0). By doing this, we immediately 
find two different regimes, which crucially depend on the 
initial value of the couplings for Hubbard-like interac¬ 
tions. To understand this, we must compare two different 
energy scales which emerge in the present problem; the 
first one is the critical energy scale Ac set by the couplings 
that diverge in the low-energy limit and the other scale 
is the low-energy cutoff of the particle-hole and particle- 
particle bubbles (say, at the wavevector q = Qx) defined 
by Eqs. © and ©, i.e., Vpkc{l—ta.n6ii)/X. Eor simplic¬ 
ity, we will consider here the renormalized Fermi surface 
that emerges at two-loop RG order to make a qualita¬ 
tive estimate as to which of the above scales win over the 
other in the two physical regimes mentioned above. 

In the first regime, which occurs for weaker couplings 
in the present model (i.e. gi < gc, where gc is a critical 
coupling parameter), it is easy to verify that Vpkc > Ac. 
This means that the low-energy cutoff in the particle- 
hole and particle-particle bubbles at the incommensurate 
wavevector q = Qx will prevent those polarization func¬ 
tions in the corresponding susceptibilities to diverge in 
the scaling limit. In this way, both PDW and GDW or¬ 
ders will remain subleading with respect to BDW and 
SSG, since these latter two orders clearly do not possess 
an infrared cutoff in their respective bubbles [see Eqs. 

and ©] for the linearized dispersion assumed in the 
present model. This result is clearly in agreement with 
several results obtained using different approaches con¬ 
cerning the spin-fermion mode l^°i^^ and it has also been 
obtained by two of us within the context of the present 
fermionic hot-spot model in a previous worb^ and, sub¬ 
sequently, confirmed by other researchers in Refj^. 

By contrast, in the second regime, and for our present 


purposes the most interesting one, the behavior of the 
system changes dramatically. This occurs for moderate 
couplings in the present model (i.e., gi > gc)-, since in 
this case it turns out that Ac > Vpkc- Most response 
functions calculated in this work will exhibit a diver¬ 
gence within such a “high-temperature” regime, which 
is due to the divergence of the couplings at one-loop 
RG order explained in the previous section. Gonse- 
quently, the PDW and GDW response vertices will now 
behave near the energy scale Ac as a power-law given by 
Am (A) ^ (A — Ac)"™, with the exponents given approx¬ 
imately by 


apDW(Q,,(„)) = -A(C3t -I- C3„) w -0.37, (19) 

«CW(Q.(.)) = ^ -0.38, (20) 

= ^(2C3v - Cst) ^ 0.38. (21) 

To establish a precise comparison among the different 
possibilities of ordering tendencies in the present model, 
we quote the previous one-loop results obtained by two of 
us in a previous worb^ regarding the critical exponents 
of AF, BDW and SSG, which are given by 

CHAF = ~-ziC2 + C2x +^3-1- C^x + ss —0.75, 

( 22 ) 

- C2c + Cap - 2C3x) ~ -0.16, (23) 

= ^(<^2c + Gic - Cix - C2x) - -0.19. (24) 

We note from the above one-loop RG results that the ex¬ 
ponents associated with GDW and PDW response func¬ 
tions (both with a d-wave form factor) at the experimen¬ 
tally observed incommensurate wavevectors Qx(y) turn 
out to be negative, which implies that the two quantities 
diverge close to the critical scale Ac at this order. The 
only exception is the GDW order at wavevectors Qx(y) 
with an s-wave form factor whose exponent is positive 
and therefore the corresponding response vertex instead 
renormalizes to zero near Ac. This latter order parameter 
turns out to be irrelevant in this case and, for this reason, 
we shall neglect this order from here on. In addition to 
this, we find, quite surprisingly, that the exponents asso¬ 
ciated with d-wave GDW and PDW at Qx(y) turn out to 
be larger than the exponents related to both BDW and 
SSG orders within this regime. We mention here that 
this qualitative trend will also hold at two-loop RG level. 
Therefore, we establish here that the PDW and the ex¬ 
perimentally relevant GDW indeed emerge as leading or¬ 
dering tendencies compared to BDW/SSC in the present 
fermionic hot-spot model. With further hole doping, the 
absolute value of the wave vectors Qx(y) that determine 
the modulation of the PDW and GDW orders will nat¬ 
urally decrease. This dependence of the observed GDW 
wave vector on doping has some support from the exper¬ 
imental situation^ and is suggestive that the Fermi sur¬ 
face is indeed playing some role in non-Lanthanum-based 



compounds. Another very interesting feature coming out 
of the above one-loop RG result is the reasonable prox¬ 
imity of the critical exponents associated with PDW and 
CDW. As will become clear shortly, this is not a coin¬ 
cidence and in fact is a precursor of the appearance of 
a new emergent symmetry at low energies, which relates 
these two orders in the present model. Finally, we point 
out that the above divergences of the response functions 
turn out to be an artifact of the one-loop approximation 
and, as we will see next, two-loop order terms will elim¬ 
inate such singularities altogether. For this reason, we 
now move on to the two-loop RG case that we are, of 
course, mostly interested in. 


B. Two loop RG 

As we have shown earlier in this paper, the renormal¬ 
ized couplings at two loops approach asymptotically only 
one IR-stable nontrivial fixed point (the HSFP), which 
controls the universal physics of the model at low ener¬ 
gies. By integrating Eqs. (Hm-lUTl) in the vicinity of this 
fixed point, we obtain that the calculated response ver¬ 
tices must obey power-laws described by (A) ~ A®"™, 
where the two-loop critical exponents are given by 

’^PDWiQ^^y-,) = -H93t + 93v) +V*^ 3.22, (25) 

= ^(- 253 . + 9*3t) 3.22, (26) 

^af = -\{9*2 + 9*2. + 9*3 + 9*3. + 253 .) + 9*^ 0.04, 

( 27 ) 

d — WQVG ^ /o * sis I Tt' O * \ I * . A OA 

^BDW — 9 (2^1c “ 92c + 9zp ~ ^ ~ 0.30, 

( 28 ) 

= 2^92c + 9ic - 9i. - 92.) + V - 0.30. (29) 

From the above result, we conclude that at two-loop RG 
level the critical exponents of all response functions cal¬ 
culated in this work become positive and, as a conse¬ 
quence, instead of a divergence displayed at one-loop or¬ 
der, they now clearly scale down to zero at two loops. 
Physically speaking, this happens because of the two- 
loop self-energy feedback (i.e., the pseudogap instability) 
onto the RG equations that acts as a “preemptive or¬ 
der” that self-consistently gaps all the correlations out, 
thus finally generating only short-range orders in the 
present model. The two-loop-order corrections also im¬ 
plies that the quasiparticle weight tends to be nulli¬ 
fied at the hot spots in the present model and that both 
uniform spin and charge susceptibilities tend to become 
suppressed in the low-energy limit, which was previously 
obtained in Ref4i. This may either point to a partial 
truncation of the Fermi surface at the hot spots (e.g., 
Fermi arcs) or it might also lead to a full reconstruction 
of the Fermi surface into pockets. Thus, the properties 
of the critical theory obtained here at two loops capture 


some aspects of the phenomenology exhibited by the un¬ 
derdoped cuprates at high temperatures. Despite this, 
we note from Eqs. (nnn-dii) that the dominant anti¬ 
ferromagnetic spin fluctuations exhibit a large but finite 
correlation length ^af- This latter fact, however, dis¬ 
agrees with the experimental situation, where the cor¬ 
relation length of AF spin fluctuations is observed to be 
only of the order of a few lattice constants in underdoped 
cuprates. 

Another very interesting property that we obtain 
straightforwardly at two-loop RG level concerns one addi¬ 
tional emergent SU (2) pseudospin symmetry that shows 
up naturally in the present model. Here we use the word 
emergent in the sense that the original fermionic hot spot 
model that we analyze in the present work is not invari¬ 
ant under this given symmetry, but only at the novel two- 
loop fixed point we obtain that this symmetry appears. 
The first emergent SU{2) symmetry that relates the d- 
wave SSC order and the d-wave BDW order parameter at 
the incommensurate wavevectors Q = (iQojiQo) was 
previously emphasized in the context of the fermionic 
hot-spot model in Ref.— and this is clearly evidenced 
by the equality of those respective critical exponents. 
This symmetry is descendant of the SU{2) pseudospin 
symmetry demonstrated by Metlitski and Sachdev^ and 
explored further by Efetov et for the case of the 
spin-fermion model with a linearized dispersion. In ad¬ 
dition to this fact, we find that the critical exponents 
associated to PDW and CDW order at the experimen¬ 
tally observed wavevectors Qx{y) (both with a d-wave 
form factor) also turn out to be numerically equal to each 
other. The underlying reason for this degeneracy is that 
these latter order parameters also turn out to be linked 
by a similar SU{2) pseudospin symmetry associated with 
particle-hole transformations in the present model (one 
for each pair of hot spots connected by (tt, tt)) that maps 
the d-wave charge order at Qx{y) onto the PDW order- 
parameter at the same wavevector. In the context of 
the spin-fermion model, this has been noted by Pepin et 
al^ (see also Refi^) by including a small breaking of 
the tetragonal symmetry (orthorhombicity) at a mean- 
field calculation and also, quite recently, by Wang et 
using a Ginzburg-Landau theory, in which case 
the CDW and PDW order parameters become compo¬ 
nents of an 5'0(4) PDW/CDW “super-vector”. There¬ 
fore, we reaffirm here this result also in the present con¬ 
text of the fermionic hot-spot model (but with no break¬ 
ing of C 4 lattice rotational symmetry) by taking into ac¬ 
count important higher-order fluctuation effects. This 
agrees with the observation that the two emergent SU (2) 
pseudospin symmetries established here manifest them¬ 
selves in the present fermionic hot spot model coexisting 
with - but not simply mediated by (see, e.g. Ref—) - 
strong short-range AF spin fluctuations. Moreover, in an 
analogous way to the previous section, such an entan¬ 
gled PDW/CDW order also turns out to be subleading 
to BDW/SSC for weak couplings in the present model 
at two-loop RG level. However, as we increase the con- 
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plings to somewhat moderate interactions, we do find a 
tendency for PDW/CDW to become leading compared 
to BDW/SSC near the HSFP obtained in the present 
work. This last result would in principle agree with the 
experimental situation, in which the leading charge or¬ 
der that emerge in non-Lanthanum-based cuprates has 
indeed a modulation determined by Q 2 ,(y).Therefore, it 
seems reasonable to speculate here that this new emer¬ 
gent SU (2) symmetry relating PDW and CDW may in¬ 
deed have profound consequences for the physics of the 
underdoped cuprates, as has been recently alluded to by 
many authors in the literature^Sr— 


V. CONCLUSIONS 

In the present work, we have investigated within a 
complete two-loop RG framework the fermionic hot spot 
model relevant to the phenomenology of the cuprates, 
which describes excitations with a linearized dispersion 
in the vicinity of eight hot spots (i.e., the points in mo¬ 
mentum space in which the AF zone boundary intersects 
a putative underlying Fermi surface of the cuprate su¬ 
perconductors at low hole doping). The present model 
can be seen as descendant of the Abanov-Chubukov spin- 
fermion model that, most importantly, includes here all 
relevant interactions between the fermions and, for this 
reason, allows one to investigate on equal footing all of its 
possible instabilities from a weak to moderate coupling 
regime. 

Here, we have explicitly studied the role played by 
two types of order in the model (CDW and PDW) in 
the vicinity of the IR-stable nontrivial fixed point ob¬ 
tained at two loops. By analyzing the CDW response 
at the experimentally relevant wavevectors (AQo: 0) 

(0, ±(5o), we were able to establish that this charge or¬ 
der is short-ranged and has a predominant d-wave form 
factor, consistent with recent STM experimental^ and 
resonant x-ray scattering^. We have also focused our at¬ 
tention on the so-called PDW order, which was recently 
proposed in the literature as a potential candidate for the 
“hidden” order to describe the pseudogap phase observed 
in underdoped cupratea^Sr— , since it may lead to a sec¬ 
ondary order parameter that breaks both time-reversal 
and parity symmetries. In this respect, we have con¬ 
firmed that the PDW order with the same modulation 
given by Qx(y) emerges as an SU (2)-degenerate counter¬ 
part of the CDW, which bears some resemblances with 
the results in the spin-fermion model obtained by Pepin 
et ali^ at a mean-field calculation and also by Wang 
et alM. In such a case, the PDW and CDW order pa¬ 
rameters should become components of an SO(4) “super¬ 
vector” (hence the denomination PDW/CDW proposed 
in Refsi^i^). 

For weaker couplings in the model, the PDW/CDW 
composite order (with a predominant d-wave form factor) 
is always subleading compared to BDW/SSC. This is in 
agreement with several works in the literature that either 


consider a spin-fermion model with weak interaction a^^i^^ 
or start from the full lattice problem and analyze it 
within a mean-field approximatio n^^! . By contrast, as 
we increase the coupling of the model towards moderate 
values, the entangled PDW/CDW order becomes leading 
compared to BDW/SSC, which clearly agrees with the 
experimental situation. We point out that our present re¬ 
sult resonates with a recent proposal by P. A. Lee^ that 
the CDW order experimentally observed in the under¬ 
doped cuprates could be interpreted in terms of a PDW, 
which could arise as a potential strong coupling insta¬ 
bility in a theory of high-Tc superconductivity from the 
perspective of doping a Mott insulator. Therefore, our 
present conclusion may also be seen in qualitative agree¬ 
ment with this point of view. 

Finally, we point out that it is crucial to implement the 
RG method for the fermionic hot model at least at two- 
loop order or beyond, since we have shown here that the 
first nonanalytic contribution to the self-energy (which is 
responsible for many effects such as, e.g., the renormal¬ 
ization of the single-particle renormalized Green’s func¬ 
tion and related quantitie o^^'^^ , the renormalization of 
the FS, and others) only emerges at two loops. Phys¬ 
ically speaking, this also happens because the two-loop 
self-energy feedback (i.e., the pseudogap instability) onto 
the RG equations acts to self-consistently gap all the cor¬ 
relations out, thus eventually generating only short-range 
orders in the present model. The two-loop-order correc¬ 
tions in addition imply that, for moderate interactions, 
the quasiparticle weight tends to be nullified at the hot 
spots in the present model and that both uniform spin 
and charge susceptibilities tend to become suppressed in 
the low-energy limit, which was previously obtained in 
Ref4i. This may either point to a partial truncation of 
the Fermi surface at the hot spots (e.g., Fermi arcs) or 
it might also lead to a full reconstruction of the Fermi 
surface into pockets. Thus, the properties of the criti¬ 
cal theory obtained here might be relevant to describe 
some aspects of the phenomenology exhibited by the un¬ 
derdoped cuprates at high temperatures. Further work 
along the lines of the one presented here would include 
to consider more realistic models in which we couple, for 
instance, all the hot spots to both the so-called “luke¬ 
warm” spots (i.e., the points that are reasonably close 
to the hot spots) and also to the cold spots in the nodal 
direction to analyze how robust are the present results 
with respect to adding more degrees of freedom in the 
system. Also, another direction that clearly deserves a 
future investigation will be to calculate transport proper¬ 
ties in the present model, such as electrical resistivity and 
thermal conductivity. Since the quasiparticles are not 
well-defined in the present model, one must use instead 
the so-called memory matrix formalisnii^ that, most im¬ 
portantly, does not assume the existence of quasiparti¬ 
cles at low energies. This calculation of non-equilibrium 
properties has been initiated in recent years to discuss 
many strongly correlated systems such as, e.g., the the¬ 
ory of an Ising-nematic transition out of a metallic state 
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which breaks the lattice rotation symmetry but preserves 
translational symmetry^, a spin-liquid model with a 
spinon Fermi surface coupled to a 17(1) gauge field^, 
and lastly the theory of spin-density-wave quantum crit¬ 
ical metals^i^. For this reason, we also plan to perform 
such an important investigation for the present fermionic 
hot-spot model in a future publication. 
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Appendix A 

In this appendix, we separate, for simplicity, the RG 
flow equations for the couplings gis, gir, gib, and (? 2 / 


that we display here up to one-loop order taking into 
account the effect of the renormalization of the FS. They 
are given by 


A-^^ — {gic + g2c)9is + {gix + g2x)gir, (Al) 

~ ^9lc + 92c)9lr + (<?lx + 92x)gis, (A2) 

= A<,;„ (A3) 

A^(2g2/ - 51b) = 0- (A4) 

As for the complete RG flow equations up to two loops 
(also taking into account the effect of the renormalization 
of the FS) of the other couplings in the present fermionic 
hot spot model described by the function Pi = Adg^/dA, 
they are given by the following expressions 


fdi = 9 i+ glx + + glp ~ 9ix92x — 93p93x — ‘^93v93t + ■:^{gix92x — ffL ~ 93p93x)9ic 

+ 2^9Ic + ffl + 9lx + 92x ~ 9lxg2x — 93p93x + glp + 93 x)91, 

/52 = 2 (ffi “ 92x — 93 ~ 93x) — ‘^glv + ^ (si + 9ic9ix + 9i9ic) + ^ (^92 - 91)93 
+ ^ [(2ff2c — 9lc) (gsp + gsx) — 2ff2cff3p53a: + 2(^2 (ffsp + “ S3p53x)] 

+ ^{92 - 92c){glx - 9ixg2x) + ■^92glx + 

/^3 = (51 — ‘^ 92)93 — 9ix{g3x — “^gsp) — 92x{g3p + 93x) + 4((73t — g3v)g3t — 2g3„ + ^[(51 ~ 2^2)^ 

+ {gic — 2g2c)^ + “^glx + ‘2glx — ‘^gixg 2 x — ‘^g3p93x + ‘^gip + ‘^glx + 93 ) 93 , 

P3t = (2gi -92 + 93 + ‘^gix - 92 x + 2g3p - g3x)93t - (51 + 53 + gix + g3p)93v + ‘^ggst, 

P3v = —{92 + 93 + g2x + g3x)g3v + ‘2rig3v, 

Pic = gic + 9 ix 92 x + gls + glr + glx - 93pg3x + ^ {gixg 2 x - glx - 93pg3x) gi 

+ 2 i9ic + 9i + 9ix + 92 x ~ 9 ixg 2 x + glp + glx — 93pg3x) gic, 

P2c = 2 (5lc + 9lx + 92 x + ‘^9ls + 25lr “ glp) + ^ (giglx + 5lc5l + 5lc) 

+ 2^92c — 52) (5L + glx ~ 9lxg2x) + ^ (252c — 5lc) 53 
+ ^ [(252 - 51) {glp + glx) ~ 25253p53a; + 2g2c {glp + glx — g3p93x)] , 

Pix = 5ic52x + g2cgix + 25is5ir + 2gixgi — 92x91 — 91x92 + (53p — 53x)53 + 4,glp — ‘^^ 93 x 931 

, V O 9% glx) 

+ 2 (5152c + 5ic52 - 2g2c52-^- Y j 


(A5) 


(A6) 


(A7) 

(A8) 

(A9) 


(AlO) 


(All) 


9 ix 


(A12) 
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FIG. 6. Relevant Feynman diagrams for the vertex corrections up to two-loop order in the present model. Solid lines de¬ 
note noninteracting fermionic single-particle Green’s function, while the wavy lines correspond to the renormalized coupling 
constants. 


P2x = gicQlx + 52c52x + 2gisgir — 3252® — 3353® “ 33® + 2 [^IcQigix — 2gicgig2x + 3lc3232® 
+ 3 i32c32® — 232c 3232® — 3i®33p33® + 32®33p33® ~ 2^2®(33® + + 2332a;, 

l^3p = (231 — 32c)33p + 31®33 + 433t “ 5233p — 32®33 - 3133® — 433„33t + - \ 2g2cg2g3p 


2 

+ 32®53® — 3i32c33p — 3ic3233p — 3i®32®33® — 3ic3i33® — 


2333p, 


1^3® — ( 231 c — 32c)33® — 3lc33p - 32®33 — 3233 ® — 33® + 2 ( 232c 3233® + 32®53p — 3 i32c33® 


3ic3233® — 3i®32®33p — 3ic3i33p 


3 L 53 ® 


2333® 


(A13) 


(A14) 


(A15) 


r 


where we have omitted the superscript i? in the renormal- up our notation and 3 is the two-loop anomalous dimen- 
ized dimensionless couplings of the model to not clutter sion contribution. In Fig. 6 , we show the corresponding 

Feynman diagrams for this calculation up to two loops. 


* hermann_freire@ufg.br 

^ T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
W. N. Hardy, R. Liang, D. A. Bonn, and M.-H. Julien, 
Nature 477, 191 (2011). 

^ T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. Berthier, 
P. L. Kuhns, A. P. Reyes, R. Liang, W. N. Hardy, D. A. 
Bonn, and M.-H. Julien, Nat. Gommun. 4, 2113 (2013). 

® G. Ghiringhelli, M. Le Tacon, M. Minola, S. Blanco- 
Ganosa, G. Mazzoli, N. B. Brookes, G. M. De Luca, A. 


Frano, D. G. Hawthorn, F. He, T. Loew, M. Moretti Sala, 
D. C. Peets, M. Salluzzo, E. Schierle, R. Sutarto, G. A. 
Sawatzky, E. Weschke, B. Keimer, and L. Braicovich, Sci¬ 
ence 337, 821 (2012). 

^ A. J. Achkar, R. Sutarto, X. Mao, F. He, A. Frano, 
S. Blanco-Canosa, M. Le Tacon, G. Ghiringhelli, L. 
Braicovich, M. Minola, M. Moretti Sala, C. Mazzoli, R. 
Liang, D. A. Bonn, W. N. Hardy, B. Keimer, G. A. 
Sawatzky, and D. G. Hawthorn, Phys. Rev. Lett. 109, 





























12 


167001 (2012). 

® J. Chang, E. Blackburn, A. T. Holmes, N. B. Christensen, 

J. Larsen, J. Mesot, R. Liang, D. A. Bonn, W. N. Hardy, 

A. Watenphul, M. v. Zimmermann, E. M. Forgan, and S. 

M. Hayden, Nat. Phys. 8, 871 (2012). 

® D. LeBoeuf, S. Kramer, W. N. Hardy, R. Liang, D. A. 

Bonn, and Cyril Proust, Nature Physics 9, 79 (2013). 

^ J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 295, 466 
( 2002 ). 

® M. Vershinin, S. Misra, S. Ono, Y. Abe, Y. Ando, and A. 

Yazdani, Science 303, 1995 (2004). 

® R. Comin, R. Sutarto, F. He, E. da Silva Neto, L. Chau- 
viere, A. Frano, R. Liang, W. N. Hardy, D. Bonn, Y. 
Yoshida, H. Eisaki, J. E. Hoffman, B. Keimer, G. A. 
Sawatzky, and A. Damascelli, Nat. Mater. 14, 796 (2015). 
“ K. Fujita, M. H. Hamidian, S. D. Edkins, C. K. Kim, Y. 
Kohsaka, M. Azuma, M. Takano, H. Takagi, H. Eisaki, S.- 
i. Uchida, A. Allais, M. J. Lawler, E.-A. Kim, S. Sachdev, 
and J. C. S. Davis, Proc. Natl. Acad. Sci. Ill, E3026 
(2014). 

R. Comin, A. Frano, M. M. Yee, Y. Yoshida, H. Eisaki, 
E. Schierle, E. Weschke, R. Sutarto, F. He, A. Soumya- 
narayanan, Y. He, M. Le Tacon, 1. S. Elfimov, J. E. Hoff¬ 
man, G. A. Sawatzky, B. Keimer, and A. Damascelli, Sci¬ 
ence 343, 390 (2014). 

E. H. da Silva Neto, P. Aynajian, A. Frano, R. Comin, E. 
Schierle, E. Weschke, A. Gyenis, J. Wen, J. Schneeloch, Z. 
Xu, S. Ono, G. Gu, M. Le Tacon, and A. Yazdani, Science 
343, 393 (2014). 

N. Doiron-Leyraud, G. Proust, D. LeBoeuf, J. Levallois, 

J. -B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 
and L. Taillefer, Nature 447, 565 (2007). 

S. E. Sebastian, N. Harrison, and G. G. Lonzarich, Reports 
on Progress in Physics 75, 102501 (2012). 

A. Abanov and A. V. Chubukov, Phys. Rev. Lett. 84, 5608 
(2000); Ar. Abanov, A. V. Chubukov, and J. Schmalian, 
Adv. Phys. 52, 119 (2003). 

M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075128 

( 2010 ). 

K. B. Efetov, H. Meier, and C. Pepin, Nature Physics 9, 
442 (2013). 

S. Sachdev and R. La Placa, Phys. Rev. Lett. Ill, 027202 
(2013). 

H. Meier, C. Pepin, M. Einenkel, and K. B. Efetov, Phys. 
Rev. B 89, 195115 (2014); M. Einenkel, H. Meier, C. Pepin, 
and K. B. Efetov, Phys. Rev. B 90, 054511 (2014). 

D. Chowdhury and S. Sachdev, Phys. Rev. B 90, 134516 
(2014). 

J. D. Sau and S. Sachdev, Phys. Rev. B 89, 075129 (2014). 
Y. Wang and A. Chubukov, Phys. Rev. B 90, 035149 
(2014). 

A. Allais, J. Bauer and S. Sachdev, Phys. Rev. B 90, 
155114 (2014). 

A. Allais, J. Bauer and S. Sachdev, Indian Journal of 
Physics 88, 905 (2014). 

A. M. Tsvelik and A. V. Chubukov, Phys. Rev. B 89, 
184515 (2014). 

S. Bulut, W. A. Atkinson, and A. P. Kampf, Phys. Rev. B 
88, 155132 (2013). 

D. Chowdhury and S. Sachdev, Phys. Rev. B 90, 245136 
(2014); D. Chowdhury and S. Sachdev, arXiv:1501.00002 
A. Thomson and S. Sachdev, Phys. Rev. B 91, 115142 
(2015). 


S. Bulut, A. P. Kampf, W. A. Atkinson, arXiv: 1503.08896 
(2015). 

“ P. A. Lee, Phys. Rev. X 4, 031017 (2014). 

D. F. Agterberg, D. S. Melchert, and M. K. Kashyap, Phys. 
Rev. B 91, 054502 (2015). 

E. Fradkin, S. A. Kivelson, and J. M. Tranquada, Reviews 
of Modern Physics 87, 457 (2015). 

Y. Wang, D. F. Agterberg, and A. Chubukov, Phys. Rev. 
B 91, 115103 (2015). 

For a comprehensive review, see P. A. Lee, N. Nagaosa, 
and X. -G. Wen, Rev. Mod. Phys. 78, 17 (2006). 

B. Fauque, Y. Sidis, V. Hinkov, S. Pailhes, G. T. Lin, 

X. Ghaud, and P. Bourges, Phys. Rev. Lett. 96, 197001 
(2006). 

L. Mangin-Thro, Y. Sidis, A. Wildes, and P. Bourges, 
arXiv:1501.04919 

C. M. Varma, Phys. Rev. B 73, 155113 (2006). 

P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 
A. 1. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 20, 
762 (1965). 

Y. Li, V. Baledent, N. Barisic, Y. Gho, B. Fauque, Y. Sidis, 

G. Yu, X. Zhao, P. Bourges, and M. Greven, Nature 455, 
372 (2008). 

J. Xia, E. Schemm, G. Deutscher, S. A. Kivelson, D. A. 
Bonn, W. N. Hardy, R. Liang, W. Siemens, G. Koster, M. 

M. Fejer, and A. Kapitulnik, Phys. Rev. Lett. 100, 127002 
(2008). 

H. Karapetyan, J. Xia, M. Hucker, G. D. Gu, J. M. Tran¬ 
quada, M. M. Fejer, and A. Kapitulnik, Phys. Rev. Lett. 
112, 047003 (2014). 

V. S. de Carvalho and H. Freire, Annals of Physics 348, 
32 (2014). 

S. Whitsitt and S. Sachdev, Phys. Rev. B 90, 104505 
(2014). 

J. Solyom, Adv. Phys. 28, 201 (1979). 

'^® M. E. Peskin and D. V. Schroeder, An Introduction to 
Quantum Field Theory (Perseus Books, Cambridge, 1995). 
V. S. de Carvalho and H. Freire, Nuclear Physics B 875, 
738 (2013). 

A. Ferraz, Phys. Rev. B 68, 075115 (2003). 

H. Freire, E. Correa, and A. Ferraz, Phys. Rev. B 71, 
165113 (2005). 

H. Freire, E. Correa, and A. Ferraz, Phys. Rev. B 78, 
125114 (2008). 

S. Blanco-Canosa, A. Frano, E. Schierle, J. Porras, T. 
Loew, M. Minola, M. Bluschke, E. Weschke, B. Keimer 
and M. Le Tacon, Phys. Rev. B 90, 054513 (2014). 

C. Pepin, V. S. de Carvalho, T. Kloss, and X. Montiel, 
Phys. Rev. B 90, 195207 (2014). 

T. Kloss, X. Montiel, and C. Pepin, Phys. Rev. B 91, 
205124 (2015). 

V. Mishra and M. R. Norman, arXiv:1502.02782 v^l. 

Y. Wang, D. F. Agterberg, and A. Chubukov, Phys. Rev. 
Lett. 114, 197001 (2015). 

®® D. Forster, Hydrodynamic Fluctuations, Broken Symme¬ 
try, and Correlation Functions (Benjamin, Cambridge, 
1975). 

S.A. Hartnoll, R. Mahajan, M. Punk, and S. Sachdev, 
Phys. Rev. B 89 (2014) 155130. 

H. Freire, Annals of Physics 349, 357 (2014). 

S. A. Hartnoll, D. M. Holman, M. A. Metlitski and S. 
Sachdev, Phys. Rev. B 84, 125115 (2011). 

®® A. A. Patel and S. Sachdev, Phys. Rev. B 90, 165146 
(2014). 


